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Space launch vehicles incorporate upper-level wind assessments to determine wind 
effects on the vehicle and for a commit to launch decision. These assessments make 
use of wind profiles measured hours prior to launch and may not represent the 
actual wind the vehicle will fly through. Uncertainty in the winds over the time 
period between the assessment and launch introduces uncertainty in assessment of 
vehicle controllability and structural integrity that must be accounted for to 
ensure launch safety. Temporal wind pairs are used in engineering development of 
allowances to mitigate uncertainty. Five sets of temporal wind pairs at various times 
(0.75, 1.5, 2, 3 and 4-hrs) at the United States Air Force Eastern Range and Western 
Range, as well as the National Aeronautics and Space Administration’s Wallops 
Flight Facility are developed for use in upper-level wind assessments on vehicle 
performance. Historical databases are compiled from balloon-based and vertically 
pointing Doppler radar wind profiler systems. Various automated and manual 
quality control procedures are used to remove unacceptable profiles. Statistical 
analyses on the resultant wind pairs from each site are performed to determine if 
the observed extreme wind changes in the sample pairs are representative of 
extreme temporal wind change. Wind change samples in the Eastern Range and 
Western Range databases characterize extreme wind change. However, the small 
sample sizes in the Wallops Flight Facility databases yield low confidence that the 
sample population characterizes extreme wind change that could occur. 


Nomenclature 

k = scale parameter of the Generalized Extreme Value estimate 

P u — probability level of a universal population 

P s = probability level of a sample population 

n = sample size 

u and v = east-west and north-south wind change components, respectively, m/s 
x = value of wind change in Generalized Extreme Value probability distribution function, m/s 

y = probability level in Generalized Extreme Value probability distribution function 

p = shape parameter of the Generalized Extreme Value estimate 

a = location parameters of the Generalized Extreme Value estimate 

I. Introduction 

Space launch vehicle commit-to-launch decisions include an assessment of the upper-level (UL) 
atmospheric wind environment to determine the vehicle’s predicted controllability and structural integrity 
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1 during ascent. These assessments occur at predefined times during the launch countdown based on the most 

2 recent measured wind data. However, the pre-launch measured winds may not represent the wind 

3 environment during the vehicle ascent. Uncertainty in the UL winds over the time period between the 

4 assessment and launch can be mitigated by engineering analysis of critical vehicle trajectory variables and 

5 development of statistical models that provide the required level of protection. Trajectory allowances as a 

6 function of time interval that reduce the maximum allowable values (“knockdowns”) of selected critical 

7 variables are made available during the countdown to ensure that the launch decision takes into account the 

8 uncertainty wind changes impose on trajectory assessments 1 " 2 . One method to quantify the distribution of 

9 UL wind change at a launch site is to compile measured temporal wind profile pairs. 3 ' 5 In order to provide 

10 robust engineering analyses the sample size must be large enough such that the sample population is 

1 1 representative of a universal population at required probability levels. 3,4 Historically, empirical statistical 

12 techniques of quantifying temporal UL wind change lacked sufficient sample sizes which required fitting 

13 the empirical UL wind change distributions to theoretical extreme value probability distribution functions 

14 in order to obtain estimates of UL wind change at extreme percentile levels. 3 Over the past decade, 

15 upgraded UL wind measurement systems have been operating at launch ranges and activities have been 

16 conducted to archive UL wind data in order for use in vehicle trajectory assessments. 6 " 7 Incorporating these 

17 archived UL wind profiles will reduce the uncertainty and increase confidence in the quantification of 

18 empirical extreme UL wind change percentile levels for use in trajectory assessments and also reduce the 

19 uncertainty in launch availability assessments. 

20 Historical databases were collected from three locations: National Aeronautics and Space 

21 Administration’s (NASA’s) Kennedy Space Center co-located on the United States Air Force’s (USAF’s) 

22 Eastern Range (ER) at the Cape Canaveral Air Force Station, USAF’s Western Range (WR) at Vandenberg 

23 Air Force Base, and NASA’s Wallops Flight Facility (WFF). These historical databases consist of a certain 

24 number of wind pairs, where two wind profile measurements spaced by the time period of interest define a 

25 pair. Sets of temporal wind pair databases for five time intervals (0.75, 1.5, 2, 3, and 4 hours) were 

26 generated from historical data at each location. 



1 

2 


II. Data Sources 

Data from multiple UL wind measurement instrumentation systems were utilized to maximize sample 


3 size in the databases for each location. Wind profiles could either be from measurements by a rawinsonde 

4 balloon, a high-resolution Jimsphere balloon or multiple vertically pointing Doppler Radar Wind Profiler 

5 (DRWP) systems. The rawinsondes are lofted by latex balloons and transmit the data back to a ground- 

6 based receiving system. Rawinsondes typically reach 30.5 km (100 kft) before the balloon bursts. Output 

7 from rawinsonde data is usually presented in pressure levels, which corresponds to unevenly spaced 

8 altitude levels. To use rawinsonde data for vehicle response assessments, data are linearly interpolated to 

9 30.5 m (100 ft) altitude levels to fill in gaps where wind data are not reported. Another balloon-based 

10 system, used only at the ER and WR for space vehicle support, makes high-resolution 30.5 m (100 ft) wind 

1 1 measurements through the use of a specially designed balloon known as a Jimsphere. 8 There are two types 

12 of high-resolution wind measurement systems that can track Jimsphere ascent. One system uses ground- 

13 based radar to track a Mylar-coated Jimsphere balloon while another system, known as Real Time (RT) 

14 Automated Meteorological Profiling System (AMPS) High Resolution Flight Element (HRFE), uses Global 

15 Positioning Satellite technology and a clear Jimsphere balloon to track ascent. 8 " 10 The AMPS system can 

16 also track ascent of a rawinsonde balloon, known as a Low Resolution Flight Element (LRFE). The 

17 Jimsphere balloon itself is more rigid than a rawinsonde, plus it contains roughness elements over the 

18 surface to disturb the air flow over the balloon and reduce the balloon’s self-induced oscillation during 

19 ascent. 8 The Jimsphere also contains a vent valve in order to maintain a constant volume and thus the 

20 balloon rises at a constant rate of 5 m/s (10 knts). However, maintaining constant volume limits the altitude 

21 range the balloon can achieve. A Jimsphere can typically reach between 16.7-18.3 km (55-60 kft). 8 

22 Vertically pointing DRWP systems are ground-based instruments that transmit and receive electronic 

23 pulses that can be converted to wind speed and direction. At a minimum, two oblique beams pointing 

24 slightly off-zenith at orthogonal azimuth angles must exist to determine horizontal wind velocities, and a 

25 vertical beam must exist to determine vertical wind velocities. The DRWP transmitted frequency and 

26 antenna size dictates the altitude range sampled and the sampling interval. The ER has a 50-MHz and five 

27 915-MHz DRWPs that, when their measurements are spliced together, can generate a wind profile from 

28 roughly 0.183-18.3 km (0.60-60 kft). 11 Data from the spliced profile are interpolated to a 15.2 m (50 ft) 
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altitude interval. 11 . Unlike balloon-based systems, the DRWP operates continuously, with the 50-MHz 

2 DRWP reporting measurements approximately every 5-mins and the 915-MHz DRWP reporting 

3 measurements approximately every 1 5-mins. Both DRWP systems produce wind profiles at vertical 

4 resolutions acceptable for launch vehicle assessments. These attributes yield orders of magnitude more 

5 profiles compared to balloon profiles available for developing temporal wind pair databases. The WR also 

6 has an equivalent set of DRWPs but due to time limitations data from both systems could not be integrated 

7 to generate vertically complete profiles. 

8 A. NASA Wallops Flight Facility 

9 Rawinsondes provide the only source of wind data at WFF, and two databases of rawinsonde profiles 

10 from WFF were obtained. The first was from the National Climatic Data Center (NCDC) Integrated Global 

11 Radiosonde Archive (IGRA) 12 database for the October 1963 through January 2000 period of record 

12 (POR). IGRA consists of rawinsonde observations from sites around the globe and contains wind and 

13 thermodynamic data at standard, surface and significant levels. The IGRA data for WFF consists of 

14 balloons released from the National Weather Service (NWS) location near WFF. The other database of 

15 rawinsondes was obtained directly from WFF, has a POR of February 2000 through January 2013, and 

16 consists of rawinsondes released at the NWS site and at WFF in support of mission operations. The IGRA 

17 database includes the rawinsonde data that was directly obtained from WFF personnel, which implies that 

18 no reason exists to include the IGRA data post December 1999 in the WFF wind pairs database generated. 

19 B. USAF Western Range 

20 Archived data from rawinsondes and Jimspheres were available for developing the WR wind pair 

21 databases. The data came from three sources: IGRA rawinsondes from January 1965 through January 2013, 

22 WR rawinsonde profiles from February 2008 through April 2012, and a Jimsphere database from January 

23 1965 through September 2001. The time overlap between the IGRA and WR rawinsonde databases was 

24 necessary because the IGRA database contains only WR rawinsonde data. Whereas, the WR rawinsonde 

25 archive contains real-time AMPS (RTAMPS) LRFE wind data in addition to rawinsonde data that is 

26 archived in the IGRA database. As stated earlier, the WR wind pair databases do not contain DRWP 

27 measurements because the DRWP systems required extensive time to process. 
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C. USAF Eastern Range 


2 Data from the ER DRWP systems provide the largest sample size and are the sole source used for 

3 database development. The spliced ER DRWP profde database has a POR of April 2000 through December 

4 2009. This POR results from the availability of quality controlled (QC’d) data for both the 50-MHz and 

5 915-MHz DRWP. 11 No rawinsonde or Jimsphere data were used because adding these data would have 

6 only increased the sample size of the ER database by 0.5%. 

7 III. Data Processing and Quality Control Procedures 

8 

9 Extensive QC of wind profde data were required to remove suspect data in individual profdes as well as 

10 in profde pairs. Automated and manual QC checks were applied on the data from each measurement 

1 1 source. The automated QC checks differed between the measurement sources and consisted of general and 

12 study-specific checks. The latter checks were necessary because the general QC checks only evaluated data 

13 in single profdes and did not check consistency within a profde pair. A profde was rejected if it failed any 

14 given QC check. The following sections present details of the QC checks for each measurement system. 

15 A. Rawinsonde 

16 Rawinsonde data from all sources went through a two-step QC process. The rawinsonde data obtained 

17 from IGRA already were processed through a set of QC checks. 12 An additional set of more stringent QC 

18 checks were performed on all rawinsonde data based on a manual review of the IGRA data. The QC checks 

19 assess data from an individual measurement perspective as well as at vertically adjacent levels. For a 

20 profde to be deemed acceptable for use, it must pass several QC checks. At least ten altitudes that contain 

21 either wind or thermodynamic data must exist. The minimum altitude in the profde must be positive. Also, 

22 all the heights reported must increase. For wind measurements the mean wind speed over the entire profde 

23 must be less than 51.4 m s" 1 (100 kt). The vector differences between wind measurements at vertically 

24 adjacent levels must be less than 51.4 m s' 1 (100 kt). A vertically adjacent check was also performed for 

25 ambient temperature where the difference magnitudes between vertically adjacent ambient temperature 

26 measurements in the lowest 3.048 km (10 kft) of the profde must be less than 20°C (68°F). Lastly the 

27 derived dew point temperature must correspond to ambient temperatures at least -60°C (-76°F). 



1 The development process first extracted wind pairs containing profiles that passed all initial QC 

2 checks. However, a rigorous, iterative process applied several more QC checks on the data to ensure 

3 sufficient data quality within profile pairs. This process removed a wind pair that contained one or two 

4 profiles that failed a QC check. Wind data with each profile in the pair must reach a minimum of ~6.1 km 

5 (20 kft). The minimum altitude requirement of 6.095 km (20 kft) was based on the minimum altitude 

6 required in order to perform certain launch vehicle assessments but insufficient for launches to altitudes of 

7 maximum aerodynamic pressure. A vertically adjacent check of wind component change was necessary 

8 since the rawinsonde data was linearly interpolated to fill in missing data that could result in artificially 

9 generated wind change. The wind change at vertically adjacent levels must not exceed 0.15 s" 1 . A related 

10 QC check was to reject profiles not containing at least 50% of the possible wind data. This check 

1 1 eliminated the potential of having an artificially large temporal wind change at the same altitude over the 

12 two profiles due to a large data interpolation in one of the profiles. A percentage check as opposed to an 

13 altitude range check was used due to the IGRA data being reported in pressure levels, which results in 

14 unequal altitude intervals and large data gaps inherently existing in valid profiles. The duplicate pair QC 

15 check was necessary for the WR wind profiles due to using multiple sources of wind profile databases with 

16 overlapping POR. If the profiles were less than or equal to 15 minutes apart they were grouped and the 

17 profile reaching the highest altitude was included in the database. Those unique profiles were then merged 

18 with the existing IGRA and RTAMPS database. 

19 The last step entailed manually inspecting each pair. The development process implemented this step 

20 after a review of the maximum wind component change and probability distributions, independent of 

21 altitude, for each pair time interval. Temporal wind change analyses have shown that wind change extremes 

22 are typically correlated to time separation, where longer time intervals imply larger extreme wind change 

23 magnitudes. 13 However, the WFF and WR wind component change at probability levels greater than 95% 

24 in the 0.75 and 1.5-hour pairs were ~50% greater than the corresponding maximum wind change at the 

25 same probability level in the 2, 3 and 4-hour pairs without manual QC. Manual inspection of the WFF and 

26 WR wind pairs for all time intervals revealed these questionable wind change values were associated with 

27 profde pairs occurring around 0000 UTC in data obtained from the IGRA database. 12 This characteristic 

28 appeared in the 2, 3 and 4-hour wind pairs for both WFF and WR. The differences observed in the 



1 questionable wind profiles seemed more characteristic of diurnal-scale wind change as opposed to short- 

2 time period wind change, which led to questioning the time stamp of each profile in the pair. An 

3 independent source provided a comparison to events of questionable wind change. Data from the National 

4 Centers for Environmental Prediction (NCEP) North American Regional Reanalysis project at the time 

5 period of interest were evaluated to determine if a large gradient in the winds existed over the time period. 

6 The majority of these cases did not corroborate with the NCEP data and, as a result, the wind pairs 

7 occurring near 0000 UTC that contained wind profiles from the IGRA source were removed from the 

8 generated database. The resultant temporal wind change distributions were better correlated as a function of 

9 time separation. 

10 B. Jimsphere 

1 1 The WR wind pairs include wind profiles from both systems that use the Jimsphere balloon. Jimsphere 

12 wind profiles were generated during launch vehicle operations and were manually QC’d by technicians 

13 prior to distribution to launch vehicle operators. The manual QC checks were performed to remove suspect 

14 data for use in flight vehicle assessments. 6 A profile was removed if its lowest altitude exceeded 122.2 m 

15 (400 ft) or if the profile contained any decreasing altitudes. All variables (altitude, wind speed, and wind 

16 direction) were removed at a specific altitude if at least one variable were missing and the remaining data 

17 above were removed after the first altitude containing missing data. Any linear wind component 

18 interpolations at the top of the profile were also removed. 

19 The selected wind pairs for the WR can be made up of two Jimspheres, two rawinsondes, or a 

20 Jimsphere and a rawinsonde. The issue with the Jimsphere/rawinsonde combination is that a difference 

21 exists in the smallest resolvable wavelengths between these two wind profiles due to their sampling 

22 intervals. The small-scale wavelengths were removed through a filtering algorithm in order to maintain an 

23 equivalent effective vertical resolution between the rawinsonde and Jimsphere systems. 8 A 244.3 m (800 ft) 

24 filter was applied to the Jimsphere based on a power spectrum analysis of the rawinsonde data shown in 

25 Fig. 1. Filtering the Jimsphere data was necessary to use wind profiles from either system interchangeably 

26 in assessing wind effects on vehicle performance. 8 The RTAMPS mean normalized power spectral density 

27 likely contains additional noise in the 152-610 m (500-2000 ft) wavelength range. The additional noise was 
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not removed from the RTAMPS data contained in the database since filtering the Jimsphere data to ~6 1 0 m 



1 (-2000 ft) would remove valid spectral content that is necessary to assess wind affects on vehicle 

2 performance. 

3 C. DRWP 

4 The ER wind pairs consist of profiles from the DRWP systems, which are designed to operate 

5 continuously with limited manual QC processing. Various algorithms and methodologies were applied in 

6 the QC process of the 50-MHz DRWP (D-50) data and in the 915-MHz DRWP (D-915) automated QC 

7 process. 7,14 Additional QC on the D-915 data included removing profiles with duplicate timestamps, filling 

8 temporal data gaps greater than 15 minutes, and checking for correct altitude progression. The resultant 

9 POR of temporal overlapping QC’d D-50 and D-915 data extends from April 2000 to December 2009. 

10 The next step after data QC entailed creating a single profile by splicing the two DRWP profiles at the 

1 1 altitude where the top of the D-915 profile, 6,100 m (20,013 ft) at most, and the bottom of the D-50 profile, 

12 approximately 2,500 m (8,202 ft), meet or overlap. Before generating the spliced DRWP profiles, the 

13 temporal and spatial (vertical) criteria are determined and applied to all DRWP profiles because the 

14 individual profiles had to match in both domains before splicing. 11 The D-50 archive contains 

15 measurements at 150 m (492 ft) intervals every five minutes prior to an instrument upgrade in 2004 and 

16 145 m (475 ft) intervals every three minutes thereafter. The D-915 archive contains measurements at 101 m 

17 (331 ft) intervals roughly every 15 minutes. The spatial criterion for the blended profiles is 50 m (164 ft) 

18 with the temporal criterion based on the time interval of the D-50 measurements. 

19 An algorithm was applied to the data to splice the D-50 and D-915 profiles into one using each of the 

20 five D-915s. First, for each D-50 timestamp, the algorithm found the closest corresponding timestamp in 

21 each of the D-915 profiles. Next, the algorithm interpolated the D-50 and D-915 profiles to 50 m (164 ft) 

22 spacing in the altitude range 0.100-18.6 km (328-61,024 ft), which corresponds to the lowest observation of 

23 the D-915s and the highest observation of the D-50. Then, the algorithm counted the number of data gaps 

24 from each instrument, flagged 300 m (984 ft) data gaps with the D-50 and 500 m (1640 ft) data gaps with 

25 the D-915, and linearly interpolated wind components within un-flagged gaps. Both profiles contained data 

26 placeholders at the same altitudes, with a transition region existing between the two profiles typically at 

27 altitudes around 2,000-3,000 m (6,562-9,843 ft). The algorithm then spliced the two profiles using a 

28 methodology that varied slightly based on the data coverage of the two profiles within the transition region. 



1 If a D-915 profile overlapped the D-50 profile, then the algorithm combined the D-50 and D-915 wind 

2 components within the transition region using a weighting scheme that provided greater weight to the 

3 D-915 at lower altitudes and greater weight to the D-50 at higher altitudes. If a D-915 profile did not reach 

4 the D-50 profile, then the algorithm linearly interpolated the winds between the highest altitude of the D- 

5 915 profile and the lowest altitude of the D-50 profile, provided the QC algorithm did not flag the gap. 

6 Each spliced profile contained winds exclusive to the D-915 below the transition region, derived winds 

7 within the transition region, and winds exclusive to the D-50 above the transition region. Splicing the D-50 

8 and individual D-915 profiles produced up to five DRWP wind profiles at a given timestamp in the archive, 

9 depending on the QC process and how many D-915s were operating. 

10 The individual spliced profiles were then combined to generate a single composite DRWP profile 

1 1 representing the wind environment at a given timestamp. The individual spliced profiles only differed 

12 below the lowest altitude of the D-50 profile. An algorithm was developed that generated composite winds 

13 at each altitude starting at 150 m (492 ft) and ending at 18,450 m (60,532 ft). The algorithm omitted data at 

14 the lowest D-915 altitude and the highest D-50 altitude due to ground clutter effects and questionable 

15 shears, respectively. The algorithm first computed a mean reference wind using the individual spliced 

16 profiles with valid winds at each altitude. Next, vector differences between the winds from each of the 

1 7 individual spliced profiles and the reference wind were used to derive weights corresponding to each of the 

1 8 profiles. Summing the product of the weights and wind components from all individual profiles produced 

19 the composite wind at each altitude. Above the lowest reporting altitude, the algorithm computed the 

20 reference wind as the mean of the reference wind described above and the composite wind at the previous 

2 1 altitude. 

22 A subsample of the DRWP archive was produced according to specified guidelines. The wind profile 

23 had to contain data at all altitudes from 250-6,096 m (820-20,000 ft) and the wind components were 

24 linearly reduced from altitudes below the lowest reporting altitude, which ranged from 150-250 m (492-820 

25 ft), to no wind at 0 m (0 ft). Profile tops extended as high as 18,450 m (60,532 ft). 
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1 IV. Wind Pair Development and Analysis 

2 

3 The criteria to select pairs, the resultant number of wind pairs at each location, statistical analyses of 

4 the sample sizes, and distributions of extreme wind changes are presented in the following sections. 

5 A. Criteria to select pairs 

6 Constraining the pair selection to the exact time spacing with the balloon-based WFF and WR profiles 

7 limits the number of pairs since balloons are released infrequently. Therefore, for each pair the time range 

8 was expanded by +/- 15 minutes to increase the wind pair sample size. For example, profiles spaced 

9 between 2.75 to 3.25 hours were treated as 3-hour pairs. For the ER, two profiles defined a pair if the 

10 desired time separation of the pair +/- two minutes separated the profiles’ timestamps. For example, a 0.75- 

11 hour (45-minute) pair has two profiles spaced anywhere from 43-47 minutes apart. The pair selection 

12 process used a two-minute window because a large number of DRWP profiles existed and at least three 

13 minutes existed between adjacent DRWP profiles. 

14 B. Sample Size 

15 Table 1 presents the number of pairs at each time interval and location. Depending on the time interval 

16 of interest, the WFF, WR, and ER databases contain 54-127 pairs, 366-548 pairs, and 260,878- 

17 297,490 pairs, respectively. The disparity in the magnitude of samples at the ER is due to the continuous 

18 operation of the DRWP versus the discrete measurements from the balloon systems used at WR and WFF. 

19 The WR’s history of supporting space launch operations that require frequent balloon releases attributes to 

20 the difference between the number of WR and WFF pairs. 

21 C. Statistical Analysis 

22 The most frequent application of wind pair databases is to apply the empirical maximum zonal (u) and 

23 meridional (v) wind change components of each profile into a persistence assessment to determine the 

24 effects of wind change over a specific time period on vehicle performance. 5 Therefore, a large sample size 

25 must exist in order to capture the largest range of maximum wind change possible. Several analyses were 

26 conducted to determine how well the sample population at each location characterized the wind change 

27 extremes. 

28 The analyses results quantify the distribution and the confidence bound (CB) in the empirical 
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maximum wind change from the various sample sizes of each pair set. Extreme wind change population 



1 distributions are usually non-Gaussian 15 , so an extreme theoretical probability function was used to fit the 

2 data. The generalized extreme value (GEV) probability distribution function (PDF) 16 " 17 provided a good fit 

3 of the extreme u- and v-component wind changes in each pair up to roughly the 99th percentile level. The 

4 GEV PDF is expressed by: 

5 y = f{x | k,[i, a) = (£) exp [l + k ^ j [l + k 5 (1) 

6 for £ A 0 and 1 + k > 0 

a 

7 where x represents each value in a distribution of wind changes, and k, u, and n denote the scale, shape, and 

8 location parameters, respectively, of the GEV estimate. Using the results from the GEV, 95% CB at various 

9 percentile levels were calculated using the Asymptotic Distribution of Percentiles (ADP) method. 18 The 

10 ADP equation is a function of the CB, sample size and percentile level of interest. The analysis uses the 

11 95% CB as a conservative approach to assess the range of extreme wind change for selected percentile 

12 levels. 

13 Representing the data annually in the distribution plots for all sites was necessary due to the limited 

14 sample populations from WR and WFF. Seasonal and monthly variability does exist in wind component 

15 distributions. 13 When using the data for engineering assessments, ideally, data should be grouped by season 

16 or month to ensure the “knockdown” value is based on a homogenous sample population. Flowever, a 

17 database subset representing a given month or season must contain a sample size large enough to 

18 characterize the environment to the desired probability level. The ER database, which contains enough 

19 samples for monthly assessments, was compared to the 2013 Cape Canaveral Air Force Station Range 

20 Reference Atmosphere (RRA). 19 The RRA contains statistics of wind variables derived from rawinsonde 

21 balloon measurements at the ER. Statistics are reported as a function of altitude for each month and 

22 annually. The comparison first computed the statistics listed in the RRA from the ER database, and then 

23 overlaid results from both sources. The annual results compared very well with some minor differences for 

24 an individual month as shown in Fig. 2. These differences are likely caused by the different sampling 

25 characteristics between the DRWP systems and the balloon-based systems that were used to derive the 


26 RRA. 



1 Annual distribution plots of the maximum change in wind component magnitude are presented in 

2 Figures 3-6 to illustrate differences in the CB as a function of sample population. For each plot, the 

3 cumulative probability, drawn from the probability density function 20 , is along the y-axis and the magnitude 

4 of the wind component’s change is along the x-axis. The sample size of the pairs is correlated to the width 

5 of uncertainty at the 95% CB for the highest percentile levels in the sample population. As the sample size 

6 decreases the width of uncertainty at the 95% CB increases. Figures 3 and 5 show the variation in the width 

7 of uncertainty at the 95% CB. In addition, a small probability density at higher percentile levels increases 

8 the width of uncertainty. Figure 3, the WFF 2-hour pairs plot, shows this attribute - where the 95% CB at 

9 probability levels above 95% in the v-component change significantly exceeds the width of the CB in the 

10 corresponding u-component change even though the sample sizes for both u- and v-changes are the same. 

1 1 The width of uncertainty in the CB for all the ER pairs is small due to the large sample size. An 

12 example of an ER pair fit to the CB is shown in Fig. 4. The deviation of the CB from the empirical 

13 distribution above the 95th percentile level is an artifact of the CB being calculated from the GEV function, 

14 which does not fit the empirical distribution well. Flowever, the poor fit is not an issue since the sample size 

15 is large enough to justify using the empirical percentiles for almost any flight vehicle assessment. 

16 For the WFF and WR samples, respectively, the 95% CB range of uncertainty increases as the sample 

17 size decreases. The WR 95% CB range of uncertainty, as shown in Table 2, at the sample size’s maximum 

18 empirical probability level was approximately 15.4 m s" 1 (30 kt) for both wind components in all the pairs 

19 with the exception for the 4-hour pairs. The range of uncertainty at the maximum empirical probability 

20 level for the 4-hour pairs is -40.1 m s" 1 (78 kt) as shown in Fig. 5. Because of the large uncertainty at the 

21 extreme empirical percentile in the 4-hour pairs, another approach was applied to quantify the confidence 

22 of the empirical wind change data. This approach uses a function that approximates the probability level of 

23 a sample population with a specified sample size to a probability level of the universal population. 3 The 

24 function makes no assumption to the form of the probability distribution function of the wind change and is 

25 defined as: 

26 P u = l + [(n-l)-£]p s " (2) 

27 where P,< is the probability that the sample contains the universal population at the sample probability P s 

28 and the sample size, n. Stated another way; a certain sample size is required to be P„ percent confident the 



1 sample contains the P s value of the universal population. Table 3 presents the confidence level of the 

2 universal population for various sample probability levels based on the sample size in each WR wind pair 

3 interval. For the 366 4-hour wind pairs, there is 88.1% confidence that the pairs contain the 99th percentile 

4 of wind change during this time interval. The confidence level exceeds 90% for the other four time periods. 

5 These samples are large enough for most vehicle performance applications; however, a low confidence 

6 exists that these samples capture wind changes at extreme (e.g., > 99th percentile) levels. 

7 The WFF samples contain the smallest number of pairs of the three locations. Due to the small sample 

8 sizes for each time period, the 95% CB for the observed wind change extremes have a large uncertainty. 

9 The range of the 95% CB for the v-component wind change range is at least 20.6 m s' 1 (40 kt) for each time 

10 period. The observed maximum 4-hour v-component wind change of 37.0 m s" 1 (74 kt) corresponds to a 

11 98th percentile level in the sample population’s distribution. Flowever, the 95% CB at the 98th percentile 

12 level ranges from ~ 40 to 89 kt (20.6 to 45.8 m s" 1 ) as shown in Fig. 6. 

13 Table 4 presents confidence levels of the universal population for various sample probabilities based 

14 on the WFF sample size. A 16.9% confidence exists that the 4-hour pairs contain the 99th percentile of all 

15 wind changes during this period. The confidence levels range from 10-36% at the 99th percentile for the 

16 other pairs. Due to the low confidence that the sample contains extreme wind changes and large uncertainty 

17 in the 95th confidence intervals at probability levels above 95%, the recommendation is to apply the 

18 extreme 4-hour wind component change for all time change intervals of interest in vehicle performance 

19 evaluations. Applying this recommendation produces more conservative results for shorter time periods, 

20 while increasing the likelihood of under-conservative results as the time period approaches 4-hours. 

21 Table 5 presents confidence levels of the universal population for various sample probabilities based 

22 on the ER sample size. The confidence level is 100% for all time periods. 

23 V. Conclusion 

24 

25 Temporal upper level wind pair databases were generated for incorporation into commit-to-launch 

26 decisions based on upper level wind assessments. Databases for five time intervals (0.75, 1.5, 2, 3, and 

27 4 hours) at the National Aeronautics and Space Administration’s Kennedy Space Center co-located on the 

28 United States Air Force’s Eastern Range at the Cape Canaveral Air Force Station, the United States Air 



1 Force’s Western Range at Vandenberg Air Force Base, and National Aeronautics and Space 

2 Administration’s Wallops Flight Facility were generated through use of historical data at each location. 

3 Multiple sources that measure upper level atmospheric winds at the requested sites were used for this study. 

4 Databases were compiled using wind profiles from rawinsonde, Jimsphere, and Doppler Radar Wind 

5 Profiler systems. Extensive quality control checks were applied on the data to remove unacceptable 

6 profiles. Statistical analyses of the resultant wind pairs from each site were performed to determine if the 

7 observed extreme wind changes in the sample pairs represent extreme temporal wind change. The resultant 

8 Eastern Range wind pair databases yielded sample sizes that characterize the extreme wind change 

9 environment and no restrictions on the usage exist. The Western Range wind pair database sample size is 

10 large enough for vehicle performance assessments up to the 99th percentile level. Including data from the 

1 1 Western Range Doppler Radar Wind Profiler systems would increase the sample size of temporal wind pair 

12 databases in an effort to characterize probability levels above the 99th percentile level. However, The 

13 Wallops Flight Facility database contains too small of a sample size to characterize extreme wind changes. 

14 This characteristic yields to low confidence existing that the observed extremes in each time period 

15 characterizes the extreme wind change environment. Therefore, for any vehicle performance applications at 

16 Wallops Flight Facility, the recommendation is to either apply the extreme 4-hour wind change values for 

17 all time change intervals of interest or incorporate analytical statistical techniques to fit the empirical 

18 distribution in order to characterize extreme wind change events. 
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Table 1 Sample size of wind pairs at each location. 
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4 Table 2 

5 

6 


Time Interval 
(hours) 

ER 

WR 

WFF 

0.75 

273,265 

435 

78 

1.5 

260,878 

401 

54 

2 

297,490 

548 

75 

3 

273,189 

508 

127 

4 

276,108 

366 

74 

TOTAL 

1,380,930 

2258 

408 


The WR 95% CB range of uncertainty for u- and v-wind component wind change at each 
pair’s maximum empirical probability level. 


Time Interval 
(hours) 

u- component wind 
change (knts) 

v-component wind 
change (knts) 

0.75 

28.47 

20.54 

1.5 

23.50 

24.17 

2 

24.24 

23.30 

3 

30.44 

34.47 

4 

78.25 

75.34 
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1 

2 Table 3 Confidence levels of the universal population for arbitrarily selected sample probability 

3 levels and the WR sample size for each wind pair time interval (Smith and Adelfang 1998). 
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Sample 

Probability 

Time Interval (Sample Size) 

0.75 hours 

1.5 hours 

2 hours 

3 hours 

4 hours 

(435) 

(401) 

(548) 

(508) 

(366) 

0.500 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.750 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.900 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.950 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.990 

0.9318 

0.9102 

0.9734 

0.9628 

0.8813 

0.995 

0.6400 

0.5960 

0.7592 

0.7215 

0.5466 

0.999 

0.0711 

0.0617 

0.1050 

0.0925 

0.0525 


Table 4 Confidence levels of the universal population for arbitrarily selected sample probability 
levels and the WFF sample size for each wind pair time interval (Smith and Adelfang 1998). 


Sample 

Probability 

Time Interval (Sample Size) 

0.75 hours 

1.5 hours 

2 hours 

3 hours 

4 hours 

(78) 

(54) 

(75) 

(127) 

(74) 

0.500 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.750 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.900 

0.9973 

0.9763 

0.9965 

0.9999 

0.9962 

0.950 

0.9065 

0.7592 

0.8944 

0.9886 

0.8900 

0.990 

0.1836 

0.1018 

0.1729 

0.3629 

0.1693 

0.995 

0.0584 

0.0301 

0.0545 

0.1332 

0.0532 

0.999 

0.0028 

0.0013 

0.0026 

0.0073 

0.0025 
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2 

3 Table 5 Confidence levels of the universal population for arbitrarily selected sample probability 

4 levels and the ER sample size for each wind pair time interval (Smith and Adelfang 1998). 
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Sample 

Probability 

Time Interval (Sample Size) 

0.75 hours 

1.5 hours 

2 hours 

3 hours 

4 hours 

(273,265) 

(260,878) 

(297,490) 

(273,189) 

(276,108) 

0.500 

1 

1 

1 

1 

1 

0.750 

1 

1 

1 

1 

1 

0.900 

1 

1 

1 

1 

1 

0.950 

1 

1 

1 

1 

1 

0.990 

1 

1 

1 

1 

1 

0.995 

1 

1 

1 

1 

1 

0.999 

1 

1 

1 

1 

1 
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Fig 1. Mean-normalized power spectral density (PSD) for the u-component (a) and v-component (b) 
WR Jimsphere and RTAMPS rawinsonde systems. 
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Fig 2. Comparison of ER DRWP profiles contained in 4-hr pairs database with the 2013 CCAFS 
RRA database statistics for February (a) monthly mean u- and v-component (solid line - CCAFS 
RRA; dashed line ER DRWP), June (b) monthly mean u- and v-component (solid line - CCAFS 
RRA; dashed line ER DRWP) and Annual (c) mean u- and v-component (solid line - CCAFS RRA; 
dashed line ER DRWP). 
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Fig 3. WFF maximum wind change from the 2-hour wind pairs with 95% CB for the u-(a) and v- 
component (b) wind changes. The magnitude of the wind component change is on the x-axis and 
probability is on the y-axis. The number (n) of pairs in the analysis is 75. 
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Fig 4. ER maximum wind change from the 4-hour wind pairs with 95% CB for the u-(a) and v- 
component (b) wind changes. The magnitude of the wind component change is on the x-axis and 
probability is on the y-axis. The number (n) of pairs in the analysis is 276,108. 
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Fig 5. WR maximum wind change from the 4-hour wind pairs with 95% CB for the u-(a) and v- 
component (b) wind changes. The magnitude of the wind component change is on the x-axis and 
probability is on the y-axis. The number (n) of pairs in the analysis is 366. 
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Fig 6. WFF maximum wind change from the 4-hour wind pairs with 95% CB for the u-(top) and v- 
component (b) wind changes. The magnitude of the wind component change is on the x-axis and 
probability is on the y-axis. The number (n) of pairs in the analysis is 74. 



